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ABSTRACT 



It is well known that Thomson scattering of Cosmic Microwave Background (CMB) photons in galaxy clusters introduces new 
anisotropics in the CMB radiation field, but however little attention is payed to the fraction of CMB photons that are scattered off 
the line of sight, causing a slight blurring of the CMB anisotropies present at the moment of scattering. In this work we study this 
blurring effect, and find that it has a non-negligible impact on estimations of the kinetic Sunyaev-Zel'dovich (kSZ) effect: it induces 
a 10% correction in 20-40% of the clusters/groups, and a 100% correction in ~ 5% of the clusters in an ideal (noiseless) experiment. 
We explore the possibility of using this blurring term to probe the CMB anisotropy field at different epochs in our Universe. In 
particular, we study the required precision in the removal of the kSZ that enables detecting the blurring term -t t ST/T in galaxy 
cluster populations placed at different redshift shells. By mapping this term in those shells, we would provide a tomographic probe for 
the growth of the Integrated Sachs Wolfe effect (ISW) during the late evolutionary stages of the Universe. We find that the required 
precision of the cluster peculiar velocity removal is of the order of 100 - 200 km s in the redshift range 0.2 - 0.8, after assuming 
that all clusters more massive than 10 14 tr 1 M Q are observable. These errors are comparable to the total expected linear line of sight 
velocity dispersion for clusters in WMAPV cosmogony, and correspond to a residual level of roughly 900 - 1800 r r //K per cluster, 
including all types of contaminants and systematics. Were this precision requirement achieved, then independent constraints on the 
intrinsic cosmological dipole would be simultaneously provided. 

Key words. (Cosmology) : cosmic microwave background, Large Scale Structure of the Universe 



1. Introduction 

The study of the scattering of Cosmic Microwave Background 
(CMB) photons in moving clouds of thermal electrons (like 
galaxy clusters) h as been a subject of a ctive r esearch sinc e 
the first works of ISunyaev & Zeldovichl d 19701 1 19721 |198Ch : 
IZeldovich& SvunvaevldT980T) . Indeed, this mechanism is one of 
the key two sources for intensity and polarization anisotropies in 
the CMB (the other one being associated to the p resence of grav- 
itatio nal fields, the so-called Sachs-Wolfe effect, (Sachs & Wolfe 
Il967l) ). The Thomson scattering changes the angular distribution 
of the CMB photons, partially erasing the original anisotropy 
pattern (due to the off scattering of photons propagating initially 
along the line of sight) and introducing new anisotropies if 
the electrons move with respect to the CMB rest frame. Due 
to its anisotropic nature, Thomson scattering also introduces 
linear polarization if the CMB shows an intensity quadrupole 
at the scattering place. The use of the polarization induced by 
Thomson scattering in galaxy clusters has been proposed as 
a probe for remote CMB quadrupoles. (ISunvaev & Zeldovich 
1980tlKamionkowski & Loeb Hl997tlSazonov & Sunvaevl ll999; 



However, this work will be devoted to the study of a blur- 
ring term, i.e., the term responsible for the smearing of the 
original intensity anisotropies of the radiation field. For an ax- 
isymmetric radiation field described by T{p)/To = Yji a iPi(M) 
(with p = cos 9 measuring deviations from the symmetry 
axis), the blurring term will affect all multipoles, whereas the 
aniso tropic scattering will slightly modify the quadrupole fe) 
(e.g. ISunvaev & Zeldovichlll980l:1zeldovich & Svunvaevll 19801: 
Suny aev & Zeldovichll981l) : 
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where a/ = i refers to the intrinsic cosmological dipole at the scat- 
tering place. However, if the scatterers show a peculiar motion 
with respect to the radiation field, an additional Doppler term 
must be added to the previous expression: 
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Chall inor et alj I2000I) . with implications for the Integrated 



Sachs- Wolfe effect (TSW. ICoorav & Baumann 2003) and for the 
characterization of the l arge scale density distribution in the ob- 
serva ble universe, (e.g.. lSeto & Sasaki 2000: Ab ramo & Xavierf 
I2007I) 



with n the unit vector along the line of sight and v the scat- 
terer's proper peculiar velocity. Since v • n/c is typically larger 
than AT /To (v ■ n/c - 7 x 10~ 4 in the linear approximation for 
WMAPV cosmology, AT/T ~ 10~ 5 as measured by COBE and 
WMAP), the blurring term should a priori be a typically few 
percent correction to the peculiar motion effect, (hereafter the 
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kinetic effect, kSZ, [Sunyaev & Zeldovich| [r980). However, this 
ratio does depend on the cluster's peculiar motion and angu- 
lar size, as we shall show later. At any rate, we show that this 
blurring term has potentially relevant cosmological significance: 
apart from the impact it has on accurate measurements of the 
kSZ effect, it can also provide an alternative measurement of the 
gas context in clusters. Likewise, if the Thomson scattering op- 
tical depth associated to the cluster is known, by using CMB 
observations (like those of WMAF-Q or Planclfl) this effect can 
be accurately predicted as an attenuation of the CMB anisotropy 
in the direction of the cluster. 

If detected on a set of cluster samples placed at different red- 
shifts, this blurring term can be used as a probe in situ of the 
CMB anisotropy field at those epochs. By looking at the varia- 
tions of this term at different redshift shells, one should be able 
to track the growth of new anisotropies arising at later times. In 
particular, it should enable performing tomography of the ISW 
effect, generated by the decay of the linear gravitational poten- 
tials at late epochs. 

This acquires particular relevance in the context of future 
surveys like the X-ray mission SPECTRUM-X/eROSITA0, 
whose prospects is to locate ~ 150,000 galaxy clusters on the 
sky, many of which will sit on top of high amplitude CMB 
intensity excursions where this effect can be measured more 
easily. Furthermore, we also show that this phenomenon should 
also provide stringent limits on the amplitude of the intrinsic 
cosmological dipole. 

This paper is organized as follows: in Section 2 we describe 
the blurring term within the context of Thomson scattering. In 
Section 3 we study its implications in the measurement of kSZ 
and remote quadrupoles at the position of galaxy clusters. In 
Section 4 we introduce the possibility of using the blurring term 
for tracking the growth of the ISW: we analyze the requirements 
in the cluster sample and in the peculiar velocity recovery. We 
observe the possibility of setting constraints on the cosmological 
dipole by using this effect in Section 5. Finally, in Section 6 we 
discuss our results and conclude. 

2. The scattering of CMB photons in electron 
clouds 

Thomson scattering modifies the angular pattern of the CMB 
intensity and polarization anisotropies. The source for new in- 
tensity anisotropies is associated to the p eculiar velocity of the 
gas cl oud with respect to the CMB frame (I Sunyaev & Zeldovichl 
1980), whereas in the case of polarization anisotropies is ex- 
clusively associated to the CMB quadrupole at the scatter- 
ing place. If the electron gas is at high temperature, then 
Compton scattering transfers energy from the electron plasma 
to the CMB photon field, distorting the CMB black body spec- 
trum a nd introducing frequency de pendent temperature fluctu- 
ations dSunvaev & Zeldovichl 1 19721 hereafter denoted as tSZ 
effect). The tSZ effect (and its relativistic corrections) has a 
definite spectral depen dence dSunvaev & Zeldovichl 1971 1 198 It 
lltoh & No zawa 2004), so hereafter we shall assume that it can 
be accurately subtracted. With this in mind, the change in the 



CMB intensity and polarizati on pattern when it crosses an elec- 
tron cloud can be written as dHu & Whitdl 19971) : 



1 WMAP's URL site: 

http : lambda . g sf s . nasa . gov/ product/map/ current/ 

2 Planck's URL site: 

http : / /www . esa . int/esaMI/Planck/index . html 

3 Spectrum-X/eROSITA's URL site: 

http : / /www . mpe . mpg . de/pro j ects . html#erosita 
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where the cl subscript refers to cluster position and tt denotes 
its optical depth. The kernels At,-j with i,j = T,+,— express 
the gen eration of new polarization /quadrupole at the cluster po- 
sition, dSazonov & Sunvaevl ll999). Note that these expressions 
are up to second order in the quan tities associated to the cluster 
(v n, Tt, Suny aev & Zeldovichll981l) . When considering higher 
order terms, new source terms for polarization and intensity in 
the direction of galaxy clusters arise, observing combinations of 
terms of the type tjv^, tjv^ for a single scattering and terms of 
the type t\ v r , r|v ( for double scattering, (y r and v t denote the 
radial and tangential components of the scatterer's peculiar ve- 
lociti e s with respect to the line of sight, ISunvaev & Zeldovichl 
119811: lltoh et al II2000I) ). We remind that, due to their peculiar 
frequency dependence, we are ignoring all terms associated to 
corrections of the tSZ and consider only those which are spec- 
trally indistinguishable from the primordial CMB anisotropies. 

Our interest in this paper will be focused in the intensity 
blurring term (-TTST(fi)/T()), which accounts for the fraction of 
photons that, initially propagating along the line of sight, were 
scattered it and never reach the observer. This term hence de- 
scribes the erasing of the CMB anisotropies at the scattering po- 
sition along the line of sight towards the electron cloud, since the 
Thomson scattering tends to isotropize the CMB angular fluctu- 
ations in that direction. The source of polarization, instead, is the 
local CMB intensity quadrupole, which is sensitive to the CMB 
at all directions in that position. In what follows, let us regard the 
galaxy cluster and group population as clouds of free electrons. 

3. Impact on kSZ estimations 

As mentioned above, the kSZ effect in clusters has the same 
spectral dependence of the intrinsic CMB anisotropies, and 
therefore extracting it requires the use of spatial frequency infor- 
mation. Its amplitude is directly proportional to the projection of 
the cluster's peculiar velocity along the line of sight. In linear 
theory, and after assuming that the initial rotational component 
of the peculiar velocity field is negligible (since it scales with the 
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Fig. 1. (a) Amplitude of the linear rms radial peculiar velocity in 
a WMAPV cosmology for a 2 x 10 14 rT 1 M Q cluster of galaxies. 
(b) Fraction of clusters on the sky above a given fractional level 
of correction on the kSZ induced by the blurring term. Dotted, 
dashed and solid lines correspond to a level of residual (instru- 
mental noise related) radial velocity of 60, 10 and km s _1 . 



inverse of the cosmological expansion scale factor), it is possible 
to relate the Fourier modes of the peculiar velocity with those of 
the matter density field (6^): 



.dDs(z)6 k 



(5) 



where H(z) is the Hubble function and Dg(z) is the linear growth 
factor for the density perturbations. After Fourier transforming 
into real space, it is easy to prove that the three spatial compo- 
nents of the peculiar velocity are uncorrelated, i.e., the tensor 
(v;(x)v/(x)) is diagonal (oc 5ij) for i = 1,2,30. A direct con- 
sequence of it is that studying the line of sight projected com- 
ponent of the real space peculiar velocity is totally equivalent to 
studying any of its three components. From eq.© it is easy to in- 



4 This is not the case in Fourier space, where different spatial com- 
ponents of V|j are correlated 
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Fig. 2. Fraction of the total CMB temperature rms that is gen- 
erated outside the cluster size. The solid line considers full sky 
CMB observations (7 m! -„ = 2), the dashed line a patch of roughly 
10°on a side (/„„■„ = 20). 



fer that, in linear theory, velocities will be Gaussian distributed, 
and each component will have a dispersion given by 
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where W(kR[M]) is the Fourier window function associated to 
a top hat filter of size given by the linear scale corresponding to 
the cluster mass M, R = [3p,„/(47r)]'^ 3 , and p m the matter comov- 
ing density. The linear matter power spectrum is given by P m (k). 
In Figure (QJ) we display the linear rms amplitude of the radial 
component of the peculiar velocity versus redshift in the stan- 
dard ACDM scenario. In an Einstein de Sitter (EdS) universe, 
its redshift scaling is oc 1 / Vl + Z, but in a ACDM universe its 
decay is even more slowly (it drops around a ~ 20 % in the red- 
shift range z e [0, 2]). Galaxy clusters and groups where the kSZ 
can be measured are biased, peculiar regions in the Universe. 
Different studies have attempted to quantify whether the pecu- 
liar velocities of these obj ects are actually biase d with respect to 
linear theory predictions. iBardeen et al.l dl986l) found that their 
velocity should be bette r approximated by that o f the peaks in 
a Gaussian distribution. ISheth & Diaferiol d2001h conducted a 
through study of peculiar velocities in simulations, in which they 
showed that also halo peculiar velocities depend upon the envi- 
ronment density, but even in the extreme cases the change in 
peculi ar ve l ocity w ould be smaller than ~ 70% of the average 
value. iPeell d2006l) found an average bias of ~ 30% of the halo 
peculiar velocity when compared to linear theory predictions. 

We have adopted an average bias of 30% in the halo peculiar 
velocities, and estimated the average correction to the kSZ that 
the blurring term gives rise to. For that, we have assumed that 
both the CMB temperature anisotropies and the kSZ fluctuations 
are Gaussian distributed, and computed the fraction of clusters 
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whose kSZ undergoes a correction above a certain level Ktsz = 
\5biurrlhsz\ = \8 C MB/[(v-n)/c+v res /c]\ due to the blurring effect: 



P{> KkSZ) 



v n 



d5' r 



P( s 'cmb) P 



v n 

c 



(7) 



In both cases, the function p(X) denoted Gaussian probability 
distribution function on the variable X. The term v res /c is in- 
cluded as a residual radial velocity term that avoids divergences: 
the solid line in Figure corresponds to v res = 60 km s _1 , the 
dashed line to v res = 10 km s _1 and the dotted line to a negligi- 



ble level of v l e 



^ km s . These residual radial velocities 



obey to the unavoidable level of noise present at the cluster's 
position when observed by experiments. The level of 60 km s _1 
corresponds roughly to a noise level of 1 jjK. Figure (fl} shows 
that more than 20% (5%) of the clusters undergo a correction on 
the kSZ that is larger than 10% (100 %). Obviously, for smaller 
radial velocities (and smaller values of v res ), the correction in- 
duced is larger. Approximately, the blurring amplitude can be 
expressed as \5T\ hhlrr ~ \6T cmb /o- C mb\ (r r /10~ 2 ) fiK. 

Once the background CMB temperature anisotropics and 
the gas content of the cluster are known, then the blurring ef- 
fect can be accurately predicted. Its amplitude will be larger 
for those regions hosting larger excursions of the CMB inten- 
sity field, and should dominate over the kSZ for those clusters 
having little radial relative motion with respect to the CMB. In 
the full sky, for WMAPV universe, the CMB intensity rms am- 
plitude is o~ cmb — 114 //K. This means that, for a survey of 
~ 150, 000 clusters (as the one to be conducted by SPECTRUM- 
X/eROSITA), around 10 clusters will show decrements of am- 
plitude larger than AIQtj fiK (4-<x), and ~ 400 clusters will be 
above the level of 350T7- fiK, (3-cr). 

When averaging over the cluster's area, it is of relevance to 
know what is the amplitude of the rms CMB fluctuations that 
is being blurred by the cluster. Figure (|2]i shows the fraction of 
this total rms that is being blurred by a cluster of a given size, 
i.e., the fraction of the total rms of a typical patch of the sky 
whose size is equal to the cluster size. This would be the effective 
CMB amplitude that is partially erased if, when pursuing the 
detection of secondary effects at the cluster's position, the CMB 
is integrated within the cluster's area. The solid line provides the 
case when all CMB sky available and /,„,„ = 2, the dashed line 
considers the case when only a limited patch is available, so that 
l m i„ = 20. In the first case, most of the anisotropy is coming from 
large angular scales, so the cluster size has little impact on the 
blurred intensity. For small surveys (/„„■„ = 20 corresponds to a 
patch size of 6 ~ 10°on a side), the rms CMB fluctuation has 
smaller amplitude and the cluster size becomes more important. 

We should also remark that the detection of this effect also 
provides an independent measurement of the gas content of the 
cluster, since the ratio of the blurred to the original CMB pattern 
is simply the cluster's Thomson optical depth tj. 

4. Implications for the ISW effect 

By measuring the intensity blurring term through the lines of 
sight towards clusters placed at a given redshift shell, then an 
estimate of the CMB temperature field at that redshift would 
be obtained. In Figure ([3^) we show the CMB TT (intensity) 
power spectrum as measured by observers placed at different 
redshifts; thick solid line corresponds to z = 0, dotted line to 
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Fig. 3. (a) CMB TT angular power spectrum as seen by ob- 
servers placed at redshifts z = 0,0.1,1 and 2 (solid, dotted, 
dashed and dot-dashed lines, respectively), (b) Same power 
spectra as in (a) after being free streamed to the present moment. 



z — 0.1, dashed line to z — 1 and dot-dashed line to z = 2. 
Since those observers are closer to the surface of last scatter- 
ing, the whole acoustic pattern shifts to larger angular scales 
(smaller multipoles). Furthermore, at redshifts larger than z ~ 1 
the contribution of the ISW is very small, and this is also vis- 
ible in the low I range. After the scattering, the angular pat- 
tern of the CMB anisotropies would be free streamed towards 
the observer, shifting the whole picture "back" to its standard 
position, (see Figure (Eb)). In Figure (0]l we display the free 
streaming of the CMB quadrupole multipole fl2,o(n) as seen by 
an observer placed at different redshifts into different multipoles 
aifl-s. The case of redshift z - 0.1 is displayed by solid cir- 
cles joined by a solid black line, z — 1 by red triangles joined 
by a dashed line, and z — 2 by green squares joined by a dot- 
dashed line. The further the remote observer is, the more power 
is aliased into high / multipoles. This streaming of CMB angu- 
lar anisotropies permits comparing the ISW pattern at different 
cosmological epochs on the angular/multipole scale. If the blur- 
ring term is observed through the line of sights corresponding to 
a population of galaxy clusters and groups situated at large red- 
shift, then it would provide a picture of the CMB pattern before 
the ISW arises. I.e., by observing the blurring term in galaxy 
clusters placed at different redshift shells it should be possible, 
a priori, to track the growth of the ISW effect with decreasing 
redshift. 
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Fig. 4. Proyection of the CMB quadrupole multipole 02,0 as seen 
by observers placed at different redshifts into different a^o mul- 
tipoles observed at present. Filled circles connected with a black 
solid line correspond to an observer placed at z = 0.1, red tri- 
angles by a red dashed line to an observer at z = 1, and green 
squares linked by dot-dashed line to z = 2. 
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Fig. 5. Maximum multipole Im to which a discrete cluster sam- 
ple on the full sky is sensitive (see equation (0). The ISW is well 
contained within Z < 30 - 40. 



In the future, several surveys in the optical and in the X- 
ray r ange (Pan-STARRSE DES 0, PAU-BAC0 (iBenitez et alJ 
2008), Spectrum-X/eROSITA) will probe the cosmological den- 



5 Pa n-STARRS' URL site: 

http : //pan-starrs. ifa.hawaii . edu/public/ 

6 DES's URL site: 

http : / /www . darkenergy survey .org/ 

7 PAU-BAO's URL site: 

http : //www . ice . csic . es/research/PAU/ 



sity field up to redshifts z ~ 0.9 — 1.1 with unprecedented sen- 
sitivity. Gravity relates the matter density distribution with the 
peculiar motion it gives rise to, so a good estimation for the kSZ 
in clusters and groups should be obtainable from the density sur- 
veys themselves, (this is indeed th e goal for cosmological recon- 
struction algorithms like ARGO. iKitaura & Enfilinl d2008l) ). The 
typical correlation length of the peculiar velocity field is ~ 40 
h _1 Mpc (comoving), which at z ~ 1 subtends around a couple 
of degrees. This means that, in the large angular scales where 
the ISW is present, we should expect to have a fairly high num- 
ber of uncorrelated estimates of the cluster peculiar velocity. Let 
us assume that we have, for a given redshift shell, a sample of 
N c i clusters on the sky. By looking only at the angular position 
of these clusters, we want to find out to what range of multi- 
poles we are sensitive. It is clear that, in order to sample a given 
multipole I, our sphere tracers must lie at distances smaller than 
~ n/l. For a set of uniformly distributed clusters on the sphere, 
we can assign an area of 4jt/N c i to each cluster, and hence an 
average inter-cluster separation of Q cc ^ 26 c = 4/ V^e/ (with 
8 C = 2/ -\fNci the radius assigned to each cluster area). Therefore, 
the maximum multipole to which our cluster sample is sensitive 
is given by 



7T I 



(8) 



Let us know assume that, as justified above, the noise is un- 
correlated from pixel to pixel (since the typical separation be- 
tween pixels corresponds to an actual distance that is larger than 
the typical correlation length of the kSZ). In this case, we may 
model the signal in our ith -pixel as 



s(hi) = sc*rfl(fii) + Sfcsz(n,). 



(9) 



We are focusing on the intrinsic contaminants whose subtrac- 
tion cannot be improved by additional observations at different 
frequencies and/or better sensitivities (such as the point source 
emission or the instrumental noise). Further, the symbol Sksz(nd 
in the equation above refers to the kSZ residual that remains 
after estimating the kSZ (and the peculiar velocity field) from 
the knowledge of the density distribution in that region of the 
Universe. We shall see that the requirements on the amplitude 
of these kSZ residuals are not too stringent, but must be inferred 
by other observations different to those in the CMB range. We 
are interested in the signal arising under the cluster's area, and 
therefore the application of a high-pass filter may be useful to 
minimize the contamination of the large scale signal. This filter- 
ing should leave some residuals, which are smaller for smaller 
clusters/groups, and lie typically at the level of 10 - 30 % the 
kSZ a mplitude (see Figure (6) in iHern andez-Monteagudo et alJ 
(120061) 1 

There will be other residuals due to the presence of radio/IR 
point sources, but, in any case, these residuals share the same 
statistical properties than the kSZ residuals, and they will be re- 
garded as the same: our goal is to set for them upper limits that 
enable tracking the growth of the ISW at late epochs. That is, 
we propose comparing the low I multipoles of the CMB at the 
high redshift cluster positions with the low I CMB multipoles 
measured from the whole sky. The difference must be due to the 
signal that arose between the clusters and the observer, i.e., the 
ISW. 

If clusters are homogeneously distributed over the sky, then it 
can be easily shown that the error in the estimation of a multipole 
a^m in the set of pixels/clusters equals 



A 2 [fl;,„ 



n 2 4-n 2 

Ql°~kSZ ~ TJ~°~kSZ' 

Jy c i 



(10) 
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with o~f sz the dispersion of the kSZ residual Sksz and Q c ; the 
solid angle assigned to each cluster. It is possible to make use of 
the 21 + 1 modes in order to provide an estimate of the power at 
multipole /. The error of such estimate is given by 



A[Q] - 



A 2 K ffl ] An o- z ksz 



21+1 



Nrl 21+1 



(11) 



and hence the S/N for the ISW detection at a given multipole 
reads 



,S\ 2 C ! , sw N cl (21 + 1) 



{-) = 

\Nh 



4n "Isz 



(12) 



The total S/N is obtained after adding this contribution from / 
2 up to / = Im given in equation dHJ: 



(d v !(d; 



(13) 



In Figure © we display the maximum / acce ssible by a clus- 
ter po pulation driven from the Sheth-Tormen dSheth & Tormenl 
1999) mass function in a WMAPV universe with <x 8 = 0.817. 
The width of the redshift shell equals Az = 0.1. Black circles ac- 
count for all clusters above 10 14 IT 1 M , whereas black squares 
and black triangles correspond for M„„„ = 2, 3 x 10 14 IT 1 M Q . We 
also display the maximum multipoles to be reached with cluster 
surveys obtained under the nominal sensitivity of the experiment 
eROSITA: big green circles correspond to a survey with a flux 
limit equal to 1.6 x 10 13 erg s _1 crrT 2 in the X-ray band [0.5, 5] 
KeV (all sky survey), whereas the big red squares correspond to 
a flux limit of 3.3 x 10 14 erg s" 1 cirT 2 in the same range (wide 
survey, see http://www.mpe.mpg.de/erosita/MDD-6.pdf 
for details). We remark that the ISW contain most of its power 
in the low multipole range, I < 20 - 30. 

In the top panel of Figure (O we display the actual angular 
density of those cluster populations at different redshifts, again 
for Az = 0. 1 . In the bottom panel, we display the required error 
in the peculiar velocity estimates (cr v = c crtsz) in order to 
obtain a residual that is 3 times below the ISW level (i.e., S/N 
= 3). Symbol coding is the same as in the previous plot. For 
the cluster sample of M mm = 10 14 IT 1 M we require errors in 
the radial peculiar velocity of the order of 100 - 200 km s _I 
in order to see the growth of the ISW at low redsfhits . These 
error requirements become more stringent when more massive 
cluster populations are used, and when higher redshifts are to 
be probed. They are at the level of the actual linear prediction 
for the kSZ in clusters (~ 220 km s -1 ), so a priori are not too 
stringent. However, we must remark that these are the upper 
limits for the total errors, and should account for not only kSZ 
residuals, but also all other types of possible contaminants and 
systematics. Note as well that the error on the kSZ recovery, cr v , 
is inversely proportional to the required S/N. 
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Fig. 6. Top: Angular number density of different cluster popu- 
lations in different redshift shells of width Az = 0.1. Bottom: 
Required errors in the kSZ recovery in order to obtain residuals 
3 times below the ISW level, (S/N = 3 in equation (fT~3T > ). 



(S/N) i= i = 1 with o-ksz ~ 100 km s" 1 and C\ S J ~ 4OO0uK) 2 , 
we obtain N c i ~ 4 x 10 3 . This is a relatively modest number 
of clusters: by including all groups and clusters present in the 
different redshift shells, the constraints on the cosmological 
dipole would improve even further and could eventually yield 
a detection. Already with 4,000 clusters, the constraint of ~ 20 
pK is between one and two orders of magnitude below the upper 
limit of the cosmological dipole that can be inferred from our 
modelling of the motion of the Local Group. 



5. Constraints on the cosmological dipole 

Let us assume now that we include the dipole in the analyses 
to be performed at the positions of the cluster distribution. 
The same requirements that allow tracing the ISW growth 
should permit, a priori, to impose constraints on the intrinsic 
cosmological dipole of the same order, i.e., at the level of 
few tens of pK. Indeed, if we use equation dT3T > to impose 



In this context, a natural question that arises is how the sub- 
traction of the measured dipole (which includes both the intrinsic 
cosmological dipole and the one induced by our peculiar mo- 
tion) affects our estimates of the intrinsic dipole in the positions 
of galaxy clusters. Let us assume that the total dipole measured 
on the whole sphere is 

r, — „intr , „\obs (\A\ 

a Um -a lm +a lm , (14) 
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built upon components of both the intrinsic cosmological and the 
observer's velocity dipoles. If, for the time being, we neglect all 
kSZ residuals in our subset of pixels containing clusters, ideally 
the total dipole measured at cluster positions would be 

«t=<(l-<^»+< s , (15) 
but in practice it will be 

<* <d - + < s + Z e '><' fl '»'' ( 16 ) 

/',m' 

where the last component accounts for the errors introduced by 
our discrete cluster (pixel) distribution, and {tj) is the average 
cluster optical depth. Since computing the dipole in our pixel 
subset is a linear operation on the map, the errors in the dipole es- 
timation will also be linear in the different multipoles conform- 
ing the map in the whole sphere. This leakage of power from 
different multipoles on the estimated dipole moment ( 1 , m) is de- 
scribed by the geometrical vector ey M i, which is only dependent 
on the angular distribution of the cluster set and can be computed 
with high precision. Note that we are ignoring t he relativistic 
effect s associated to the motion of the observer, (Chlub aet alj 
120051) . If, as usually done in CMB analysis, the total dipole mea- 
sured in the entire sphere is subtracted, then, in the cluster pixel 
set, we would be left with 

£<n^~< T r><£+ J] **a,„. (17) 

l'>l,m' 

I.e, the rms error on the limits for the intrinsic cosmological 
dipole will be of the order 

( Af JntrA 1/2 H/'>l,m' g/',m'Qf',m' 

l A KJJ ^rj ■ ( 18 ) 

According to this expression, it becomes critical to choose a con- 
venient cluster/pixel set that fulfills e/'>i, m '/(Tr) ^ 1, so that 
interesting limits on a'"' r can be set. If the kSZ removal is accu- 
rate and/or the large angle power of the kSZ is negligible, then 
the residuals of our IS W estimation will be dominated by those 
given in Figure 6 of IHernandez- Monteag udo et al.l d2006l) . In 
that case, for a population of ~ 500 galaxy clusters, the average 
residual dipole should be of the order of 30/ ^JNrfX V3 ~ 2 pK if 
we assume an angular size of 20 - 30 arcmins, but lower (~ 0.4 
pK) for smaller sized (~ 3 arcmin) clusters. Thes e figures ignore 
the im pact of instrumental noise. However, Kashlinskv et alj 
(2009) were able to achieve a sensitivity level of a few pK in 
each component of the dipole computed on a ~ 600 cluster set 
with WMAP data. A residual dipole amplitude of ~ 2 pK on a 
cluster set of (tt) ~ 5 x 10~ 3 would already correspond to a con- 
straint on the primordial dipole at the level of ~ 450 pK, which 
is well below the observed dipole. For the same cluster set, a 
residual amplitude of ~ 0.4 pK (as quoted for ~ 3 arcmin clus- 
ters) should yield a constraint at ~ 77 pK, which is close to the 
expected value of the intrinsic dipole (~ 30 pK). We conclude 
that ongoing CMB experiments like Planck, ACT or SPT should 
be able to set strong constraints on the intrinsic CMB dipole. 

6. Discussion and conclusions 

Current CMB high angular reso lution su rveys like ACT 
dFowler & ACT Collaboration! 120061) or SPT dRuhl et al.ll2004T) 
are scanning the sky at frequencies in which the identification of 
galaxy clusters and groups should be possible via the tSZ and 



the kSZ effects. Nominally, these experiments should be able to 
detect via the tSZ all clusters more massive than 2 x 10 14 IT 1 
M at large significance, but many smaller cluster and groups 
should remain close to the detection threshold. Once clusters 
have been identified by their tSZ distortion, attempts to detect 
the kSZ in a subset of those systems can be conducted. The 
kSZ is our only probe for peculiar velocities in the high redshift 
universe, and those can themselv es be used as probes for Dark 
Energy and missing barvons (jHernandez-Monteagudo et all 
2006; iHernandez-Monteagudo & Hoi 120091) . In this context, a 
precise characterization and correction of all possible contam- 
inants is critical. The impact of the blurring effect on the kSZ 
amounts to ~ 20% for 10% of the clusters and groups, and be- 
comes more important for those objects with small radial projec- 
tion in their peculiar motions. Provided that the Thomson optical 
depth of the cluster is known, then this effect should be accu- 
rately predicted from background CMB observations (as those 
available from, e.g., WMAP or Planck). The measurement of 
this effect is itself a measurement of the gas content of the object 
under study. 

This effect mirrors the CMB intensity at the epoch of scatter- 
ing, and this is relevant for secondary anisotropies which arise 
at late times (like the ISW effect): a detection of the blurring 
effect in objects placed at different redshift shells would pro- 
vide the picture of the growth of the ISW at recent cosmological 
times. The measurement we are proposing here is statistical, and 
therefore it does no t require a high S/N in each cluster (just in 
the same way as in IHern andez-Monteagudo & Sunvaevl {2008) 
for the kSZ - E polarization mode cross correlation). Our ar- 
guments here are therefore similar to those given in that work: 
cluster and group positions can be inferred from observations 
in a wide range of frequencies (optical, IR, X-ray, millimeter), 
many of which are using those objects for studying the nature of 
Dark Energy. The critical issue is the nature and the amplitude 
of residuals in the kSZ/CMB estimation at the cluster positions. 
To what extent do errors in the IR/radio point source subtraction 
and/or in the characterization of the local peculiar velocity field 
actually endanger the ISW measurements? This should critically 
depend on those errors being systematic or not. If those residuals 
can be regarded as independent from cluster to cluster, then the 
viability of this project should hinge exclusively on the actual 
precision to which the point source and kSZ residuals can be re- 
moved from the clusters' area. Since the ISW contains most of 
its power at low multipoles {Im < 10 - 20), targets may be chosen 
at a convenient distance in order to minimize the required sky 
coverage. A possible strategy would then consist on uniformly 
distributing N pa , c h patches on the sky (with N patc h ~ {4/n Im) 2 ), 
lying a distance 9 ~ 4/ ^N pat ch away, and scanning deep through 
each of those patches until finding a set of sources of high S/N 
at different redshifts. This would improve the efficiency of the 
survey (since a minimum amount of area would be scanned), 
at the expense however of improving the flux/mass thresholds 
shown in Figure (O. Hence one would encounter here a trade off 
between flux sensitive and sky coverage. Whatever approach is 
finally chosen, it should also provide strong constraints on the 
cosmological dipole, which would be independent from those 
imposed from the local dipole and the local velocity flows. 

In this work we propose using, by first time, the so-called 
blurring term in Thomson scattering for cosmological purposes. 
The small fraction of scattered off CMB photons which are de- 
viated when crossing a galaxy cluster/group should provide in- 
formation on how the CMB anisotropy field was at the time of 
the scattering. If those objects are far enough, then the CMB 
at that epoch should lack the ISW component that has arisen 
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recently, and this difference could a priori be picked up after 
removing all other signals present in the cluster with enough ac- 
curacy. Assuming that tSZ residuals are negligible, we find that 
a typical error of 100 - 200 km s _1 in the cluster peculiar ve- 
locity reconstruction is required for all clusters more massive 
than 10 14 h -1 M in order to trace the growth of the ISW at late 
times. These errors are comparable with the linear expectations 
for the peculiar motions of those objects, which involves that ( i) 
the blurring correction to the kSZ is general of relevance for the 
estimation of the latter and ( ii) not very precise corrections for 
the kSZ are required. Those amplitudes roughly correspond to an 
error of 900 - 1800 t t pK per cluster. The same level of errors 
would provide stringent constraints of the intrinsic cosmologi- 
cal dipole. Current and future large scale structure surveys like 
eROSITA, Pan-STARRS, DES, PAU-BAO, ACT or SPT should 
already provide enough group and cluster candidates at the rele- 
vant redshift ranges. Therefore, the critical point is the feasibility 
of an accurate enough kSZ/tSZ/point source subtraction in future 
high resolution CMB observations. 
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